Dynamical evolution of quantum oscillators towards equilibrium 
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A pure quantum state of large number A'^ of oscillators, interacting via harmonic coupling, evolves 
such that any small subsystem n << A'' of the global state approaches equilibrium. This provides 
a novel example where stationarity emerges as a natural phenomena under quantum dynamics 
alone, with no necessity to bring in any additional statistical postulates. Mixedness of equilibrated 
subsystems consisting of l,2,...,n << N clearly indicates that small subsystems are entangled 
with the rest of the state i.e., the bath. Every single mode oscillator is found to relax in a mixed 
density matrix of the Boltzmann canonical form. In two oscillator stationary subsystems, intra- 
entanglement within the 'system' oscillators is found to exist when the magnitude of the squeezing 
parameter of the bath is comparable in magnitude with that of the coupling strength. 
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I. INTRODUCTION 

Deducing the statistical distribution in many particle 
systems as an intrinsic property, resulting solely from 
quantum dynamics, has attracted much attention in the 
literature In this context, considerable interest has 
been evoked recently [1, [3, 01 on a deeper understand- 
ing of the basic mechanism of equilibration [s*] , occurring 
as a natural consequence of quantum dynamical evolu- 
tion - without invoking any additional statistical assump- 
tions. More specifically, equilibration is realized entirely 
in a quantum mechanical setting - with the key element 
ascribed to quantum entanglement between the system 
and the environment. Erstwhile statistical postulates on 
ensemble averaging over initial distributions are not re- 
quired at all - as quantum dynamics of individual pure 
states of a many body physical system itself leads to equi- 
libration of smaller subsystems. This features an excit- 
ing foundational development, where 'subjective' lack of 
knowledge in terms of statistical ensemble averaging is 
replaced by the 'objective' randomness due to entangle- 
ment f6j. 

Based on powerful general arguments, Linden et. al. 
established [j] that an overwhelming majority of pure 
quantum states of interacting large quantum systems 
evolve such that any small subsystem approaches a sta- 
tionary state. This brings out an important implica- 
tion: dynamics of almost every pure many body quantum 
state, envisages stationarity of any small subsystem as an 
inherent property, with all the statistical ingredients al- 
ready built within the basic quantum framework itself. 

In this paper, we present an explicit analysis of 
the quantum evolution of an initially uncoupled pure 
squeezed state of a large number N of oscillators, sub- 
jected to a harmonic interaction Hamiltonian, result- 
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ing eventually in equilibrium of any small subsystem of 
n << N oscillators of the global pure state - with the rest 
acting as the bath. Note that in Ref. [3| , the smallness of 
the system in relation to the size of the bath is described 
by their respective dimensions. In the present case, in- 
dividual systems constituting the whole state are infi- 
nite dimensional (being harmonic oscillators) and so, the 
smallness of the subsystems of the global quantum sys- 
tem consisting of N oscillators is expressed legitimately 
in terms of the number n << A^ of a subset of oscilla- 
tors under consideration. Essentially, we find that every 
small subsystem of oscillators tend to relax in a station- 
ary state, specified by a mixed density matrix - which is 
indeed a signature of quantum entanglement of the 'sys- 
tem' and the 'bath', emerging due to quantum dynamics. 
It may be worth pointing out here that due to its mathe- 
matical transparency and simplicity, the physical model 
of a linear assembly of coupled oscillators has played a 
paradigmatic role in understanding difficult formal prin- 
ciples underlying statistical mechanics A detailed 
investigation of quantum dynamics of pure states lead- 
ing to equilibrated small subsystems in this model would 
therefore be illuminating. 

The paper is organized as follows: In Sec. II, we de- 
scribe the physical model of A^ harmonically coupled os- 
cillators and discuss its exact solution using symplectic 
transformations - which offers a most natural elegant ap- 
proach to the problem of interest. We then identify the 
N X N symplectic transformation corresponding to uni- 
tary time evolution in this model. This is followed by 
Sec. Ill, where we analyze the time evolution of a global 
pure uncoupled squeezed state (which is not an eigenstate 
of the Hamiltonian) of the system-bath oscillators. Under 
the assumptions of continuum limit N oo and weak 
coupling approximation wc show that any small subsys- 
tem of the whole pure state of N oscillators exhibits a 
stationary long time behavior. Every single mode oscil- 
lator system is shown to relax in a mixed state of Boltz- 
mann canonical form, with effective temperature related 
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to the squeezing parameter of the bath oscillators.lt is 
also shown that any two oscillator system approach a sta- 
tionary mixed state with an intra-entanglement surviving 
whenever the magnitude of the bath squeezing parame- 
ter is comparable in magnitude with that of the coupling 
strength. Sec. IV has concluding remarks. 



II. THE PHYSICAL MODEL 

We consider a long chain of coupled harmonic oscilla- 
tors, the Hamiltonian of which is given by, 
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Here qi , pi denote position, momentum operators of the 
oscillators, satisfying the canonical commutation rela- 
tions [qi,pj\ — ihSi^j. It is convenient to define a 2N- 

component operator column ^ of dimensionless variables. 



and express the commutation relations compactly as, 

[LAp] = jr„^; a,/3= 1,2,...,27V, (2) 



where, T — 
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/ denotes the N x N unit matrix. 



A general real homogeneous linear transformation on ^ 
preserving the canonical commutation relations ([2]) is a 
27V X 2N symplectic transformation Q 5 e Sp(2A^, R) 
and there exists a corresponding unitary operator U{S) 
on the Hilbert space on which the operators ^ act: 
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Here, we restrict to Gaussian quantum states of the 
system and the bath. These are completely character- 
ized by the first and second moments of ^, arranged con- 
veniently in the form of 27V x 2iV covariance matrix V 

as, 

Va^p = i ({Ala, Aip}), a, /? = 1, 2, . . . , 2iV, 

where = i - (|), {61,62} = O1O2 + 6261 and 
(O) — Tr[/50] denotes the expectation value of the 
operator O in the quantum state p. Under symplectic 
transformation, a Gaussian state is mapped to another 
Gaussian state characterized by the covariance matrix 
V' = SVS^. 



Time evolution of the elements of the variance matrix 
under U{t) = exp{— itiJ/fi,} may be identified as a sym- 
plectic transformation (as the Hamiltonian of Eq. ^ is 
a quadratic in the canonical operators): 

V^p{t) = iTr[p(t){Ae„(0),A|^(0)}] 

= ^TT[m^Ht){Ma{0),Aip{0)}U{t)] 

or V{t) = S{t)V{0)S^{t), where S{t) denotes the 
2iV X 2iV symplectic matrix corresponding to the uni- 
tary time evolution on the Hilbert space of the quantum 
state. The explicit structure of the symplectic transfor- 
mation matrix S{t) associated with the dynamical evo- 
lution U(t) in the present model is readily identified, as 
will be outlined in the following. 

We first express the Hamiltonian ^ in the following 
quadratic form 
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where the elements of the N x N bolck matrix A are 
given by. 



A, 



(2e + 1) Sij - e [Sij+i + S^+ij] , 
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Identifying the real orthogonal transformation cr which 
diagonalizes the real symmetric matrix A i.e.. 



ctAct^ = A = diag(A((/)i), A(02), 
= 1 + 2e(l - cos (j^i - 



,K<I>n)), 
In 
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(5) 
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sin(s0i), 



we express the Hamiltonian ^ in its decoupled struc- 
ture: 



H = ^ 
2 



A2 

N 
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Pi + Q 



U{S), (6) 



where 5 = A 4 cr© A 4 is a symplectic transformation Q 
on the 2A'^ component operator column ^. 

Thus, we obtain the 2N x 2N symplectic matrix S{t) 
corresponding to the unitary time evolution operator 



S{t) 



cos{ujtA-^) A 2sm{ujtA-2) 
-A2 sm{ujt A 2 ) cos{ujt A 2 ) 



(7) 



We proceed now to investigate the quantum evolution of 
a pure uncoupled squeezed state of oscillators. 
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III. TIME EVOLUTION OF PURE 
UNCOUPLED SQUEEZED STATE 

First, we decompose the global quantum state of A''- 
oscillators into two parts: n << N 'system' oscillators 
i.e., a n oscillator subsystem and the rest of the whole 
state, the bath. We consider an initial state of the whole 
system to be pure product states of individual oscillators, 

|*(77, fi, t = 0)) = \Mv, t = 0)) \Mv, t = 0)) ® . . . 

^\Xs{fi,t ^ 0)) ® . . . ^ \Xsifi,t = 0)) 

n system oscillators 

^\Mv, i = 0)} ® . . . ® iMv, t = 0)), (8) 

where the oscillators in the bath are in the squeezed state 

|0,(ry,t = O)) = J7(5(77))|O). (9) 

Here, U{{S{r])) denotes the squeezing operator (with 
associated 2x2 symplectic matrix given by [H, 
S{r]) = diag(e~''/2, e''/^}. The initial state of each 



of the system oscillators 



\xs{pA^o))^uiSi^lm, 



(10) 



is characterized by the squeezing parameter fi. (Here, |0) 
denotes the ground state of the oscillator.) 

It may be noted that initially, the whole system-bath 
state is a product state of oscillators and the subsystems 
are not already in a stationary state, when the couplings 
are switched on at t = 0+. In other words, the global ini- 
tial pure state ^ is not an energy eigenstate of the total 
Hamiltonian - this being a trivial case leaving the sub- 
systems stationary under dynamical evolution. Also, in 
contrast to the case where the bath is initially in thermal 
state characterized by a temperature T, here the bath is 
in a pure state specified by a squeezing parameter rj. 

The initial variance matrix of the system-bath pure 
state is given by. 



{DQ{r^,ti)(BDp{r^,li)) (11) 



agonal matrices, 

DQ{■n,^J) = diag(e~'',. 



(12) 



It is convenient to split the variance matrix V^i], ii; 0) 



as, 



i=l 



where, ^(ry; 0) = - (e / © e'' /) 

Y,V^{V,^^■,0) - V{7j,^i;0)~v{Tj-:0)- (13) 

The non-zero elements of Vi{r], /i; 0) are readily identified 
as (see Eqs. (III1),(II2]), 

K(ry,M;0)]„+.^,„+., = ^(6^-6") 

Temporal evolution of the quantum state ^ is entirely 
determined by the symplectic transformation S{t) (given 
by Eq. ([7])) on the variance matrix, 

n 

1=1 



(14) 



where the blocks Dq^q, /i), and Dp{ri, fi) are x di- 



where VQQ{r], fi;t),Vpp{r], ^j,;t) and VQp{r], ^j,;t) respec- 
tively denote the N x N diagonal and off-diagonal blocks 
of the variance matrix. 

In an infinitely long chain {N —> oo), closed form an- 
alytical expressions are obtained for the elements of the 
variance matrix V{ri, fi; t) (by replacing the discrete vari- 
able '(1)1 ~ of ([5]) by a continuous parameter ' and 

the sum "j^Y^iLi ' by the integral /J" d(/) ' ): 



[yQp{r],^l]t)]^^ 



e-''cff{t) + e'^sff'\t) + 



^1 



n 



i=l 



i=l 



1=1 



'ri 



j=i 



(15) 
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where we have denoted, 



C. 



[a.K) 



sin(s0) sin(Z0) A''((/)) cos'^MA^ ((/))] 
sin(s0) sin(?(/)) A'"((/.) sin''[wiA^((?!))] 
I 



(16) 



K+2k 



These re- 



with A(0) = £7"-^ (1 + 27 cos 0), 7 
suits are formally exact in the long chain limit. To make 
their meaning evident, one resorts to the weak coupling 
approximation e « 7 << 1, in which case the following 
standard form 



)cos(s0) cos[a;(l — 7COS (/))] = Js{jx) cos(x — —) 



(where Js{x) denotes Bessel function of integral order) 
can be employed to simplify C^j'^'ht), S^'^i'^^t) of ^ 
up to order 0(7) - thus reducing the elements of the 
variance matrix T^(?7, /i; t) (given by Eq. (fT5|) ) to time de- 
pendent cos(rit) or sin(rit) functions, oscillating rapidly 
with Bessel functions Jg (7^^^) being the amplitudes (with 

n 



^tn^ y ^^^^S time behavior of the system gets 
specified by the asymptotic decay of Bessel functions i.e., 
lim Js{x) — > \x\~^ for x » s. More specifically, we 
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find. 



C, 



(2,0) 



,(1,±|) ^(1,0) 



■'s.l 



s,l 



^(2,±1) 



^{SsaTi[Ss,i+i+Ssj-i]), (17) 



in the limit jQt — > 00. 

It is thus evident that a subsystem of n << iV oscil- 
lator relaxes in a steady state, specified by the 2n x 2n 
variance matrix, V^""^ = Vg"' ® with elements. 
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1 / e~''7 
= - coshT^H —[Ss,i+i + Ss.i-i 

si Z \ Z 



1 ( e^^ 

- ( cosh?7 - -y[Ss,i+i + Ss,i-i] ] (18) 



(which exhibit a correlation {QjQj±i) = ^ and an 

anticorrelation {PjPj±i) = —-^eP between neighbors). 
Evidently, the equilibrium state of the 'system' is inde- 
pendent of its initial form (i.e., it does not contain the 
squeezing parameter ji of the initial system oscillators), 
but depends on the squeezing parameter -q of the bath. 

In particular, the variance matrix V^^^ (t) of any single 
mode subsystem of the dynamically evolving global quan- 
tum state eventually converges, in the limit t >> 
to an 'equilibrium' structure. 



2 



cosh f] 
cosh rj 



(19) 



which corresponds to a mixed density matrix of the fa- 
miliar Boltzmann form, 



^(1) 



Pr 

2m 



K 



(20) 



The inverse temperature (3 is related to the squeezing 
parameter ri of the bath via, P ~ coth^^ [cosh77]. 
The purity [ij] of the single oscillator equilibrium state 
1/(1) = Tr[(p(i))2] = [2Vdet(F(i))]-i = [coshr;]"! < 1 
captures the system-bath entanglement. This also re- 
flects in the increase of the von Neumann entropy [T^ 
of the state from its initial value zero to the equilibrium 
value 



5(p(i)) = -Tr[p(ihnp(i)] 



1-^(1) 
2i.(i) 



In 



iW^ 

1-1.(1) 



In 



2vW 



1 



/(I) 



Any two oscillator subsystem is found to eventually 
relax in a 'stationary' state specified by the two mode 
variance matrix. 



/ cosh 



y(2) = t 



\ 



2 

cosh 77 





cosh rj — 

— ^ cosh 77 J 



(21) 



The stationary density matrix of the two oscillator sys- 
tem has its purity, i/^^) = [4-\/det(F(2))]~-^ « [coshry]"^ 
which is clearly less than 1, and reveals the system-bath 
entanglement. One finds internal entanglement between 



)(cosh 77 



the two oscillators if [13| (cosh 77 

2rj 2 

^—^) — cosh2?7 < 0, and this happens when the squeez- 
ing parmeter 77 of the bath is comparable in magnitude 
with the coupling 7. Thus, one finds a tradeoff of entan- 
glement within the system oscillators and that between 
the system-bath - as survival of internal entanglement 
implies nearly vanishing mixedness of the two oscillator 
system (i.e., i'^^^ ~ 1 for small squeezing parameter). 



IV. CONCLUSIONS 

In summary, we have shown, through an explicit anal- 
ysis under the assumptions of continuum limit N 00 
and weak coupling approximation e « 7 << 1, that any 
small subsystem of the whole pure squeezed state of N 
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oscillators evolving under hamonically coupled Hamilto- 
nian, approaches equilibrium. This provides an excellent 
example in which stationary behaviour of any small sub- 
system of a large global pure state is shown to emerge 
as a natural phenomena consequent to quantum dynam- 
ical evolution - without the aid of any additional statis- 
tical postulates [4|]. Equilibrated subsystems consisting 
of 1,2, ...,n << N oscillators are found to be mixed, 
revealing their quantum entanglement with the rest of 
the system i.e., the bath. A single mode oscillator 'sys- 
tem' is shown to relax in a mixed density matrix in the 
Boltzmann canonical form. The connection between the 
squeezing parameter of the bath and the temperature is 
a new feature of our work. We also find (in the case 
of n = 2) that entanglement within the 'system' oscilla- 
tors in the steady state survives only when the squeezing 
parameter and the coupling strength are of comparable 
magnitude. These features on two and more subsystem 
oscillators, to the best of our knowledge, have not been 



recorded in the literature. 

The long chain limit and the weak coupling approxi- 
mation made our analysis amenable to analytical results. 
One has to resort to numerical approach to evaluate the 
integrals in the continuum, strong coupling limits. 
Also, the finite N limit may be addressed with the help 
of numerical investigations both in the strong and weak 
coupling limits. These issues would be of interest from 
a foundational point of view and we plan to address this 
issue in a separate communication. 
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